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Co-electrolysis  of  C02  and  H20  in  a  solid  oxide  electrolyzer  cell  (SOEC)  offers  a  promising  way  for  syngas 
production.  In  this  study,  an  electrochemical  model  is  developed  to  simulate  the  performance  of  an 
SOEC  used  for  C02/H20  co-electrolysis,  considering  the  reversible  water  gas  shift  reaction  (WGSR)  in 
the  cathode.  The  dusty  gas  model  (DGM)  is  used  to  characterize  the  multi-component  mass  transport  in 
the  electrodes.  The  modeling  results  are  compared  with  experimental  data  from  the  literature  and  good 
agreement  is  observed.  Parametric  simulations  are  performed  to  analyze  the  distributions  of  WGSR  and 
gas  composition  in  the  electrode.  A  new  method  is  proposed  to  quantify  the  contribution  of  WGSR  to  CO 
production  by  comparing  the  CO  fluxes  at  the  cathode-electrolyte  interface  and  at  the  cathode  surface. 
It  is  found  that  the  reversible  WGSR  could  contribute  to  CO  production  at  a  low  operating  potential  but 
consume  CO  at  a  high  operating  potential  at  an  operating  temperature  of  1 073  K  and  inlet  gas  composition 
(molar  fraction)  of  H20:  49.7%,  C02:  25%,  H2:  25%,  CO:  0.3%.  In  addition,  the  contribution  of  WGSR  to  CO 
production  also  depends  on  the  operating  temperature  and  inlet  gas  composition. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Hydrogen  can  be  produced  in  a  sustainable  manner,  such  as  by 
photocatalytic  water  splitting  [1  ],  gasification  of  biomass  [2],  solar 
thermochemical  water-splitting  [3],  or  water  electrolysis  driven  by 
solar  cells  or  wind  turbines  [4].  Among  all  the  above  technologies, 
water  electrolysis  is  a  practical  and  efficient  method  for  large-scale 
hydrogen  production.  Alkaline  electrolyzers  and  proton  exchange 
membrane  (PEM)  electrolyzers  usually  work  at  room  temperature. 
Solid  oxide  electrolyzer  cells  (SOECs)  use  the  same  materials  with 
solid  oxide  fuel  cells  (SOFCs)  but  work  in  a  reversed  mode  at  a 
high  temperature  (i.e.  1073K).  Compared  with  alkaline  and  PEM 
electrolyzers,  SOECs  consume  less  electricity  as  part  of  the  energy 
needed  for  water  splitting  is  in  the  form  of  heat  [5].  Because  of  their 
great  potential,  SOECs  have  received  increasing  interest  in  recent 
years  [6-20].  Various  materials  have  been  developed  to  fabricate 
SOEC  for  hydrogen  production  by  steam  electrolysis  [21-23].  Sev¬ 
eral  mathematical  models  have  been  developed  to  predict  the  SOEC 
performance  at  various  levels  [24-32]. 

In  addition  to  steam  electrolysis,  SOECs  can  be  used  to  elec¬ 
trolyze  C02  for  CO  and  02  production  [33-39].  Co-electrolysis 
of  C02  and  H20  has  also  been  demonstrated  to  be  feasible  for 
simultaneous  production  of  H2  and  CO  [33-44],  which  can  be 
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subsequently  processed  for  synthetic  fuel  production.  In  an  SOEC 
used  for  H20/C02  co-electrolysis,  3  reactions  take  place  simulta¬ 
neously,  namely  H20  electrolysis,  C02  electrolysis,  and  reversible 
water  gas  shift  reaction  (WGSR).  A  common  understanding  on 
co-electrolysis  is  that  the  reversible  WGSR  should  always  con¬ 
tribute  to  CO  production  [38,40].  However,  it  is  still  not  clear  to 
which  degree  the  reversible  WGSR  is  responsible  for  CO  produc¬ 
tion  in  the  SOEC  [45].  In  addition,  the  existing  studies  on  H20/C02 
co-electrolysis  are  all  experimental  in  nature,  with  aims  to  demon¬ 
strate  the  feasibility  of  this  technology  or  to  develop  new  materials 
for  performance  improvement.  The  present  literature  is  lacking 
detailed  mathematical  modeling  of  the  SOEC  used  for  H20/C02 
co-electrolysis. 

In  this  study,  an  isothermal  electrochemical  model  is  devel¬ 
oped  to  characterize  the  performance  of  an  SOEC  used  for  C02/H20 
co-electrolysis.  It  is  an  extension  of  the  previous  models  for  C02 
electrolysis  and  H20  electrolysis,  respectively  [29,46].  Both  the  co¬ 
electrolysis  and  the  reversible  WGSR  inside  the  SOEC  cathode  are 
considered.  The  model  is  validated  by  comparing  the  simulation 
results  with  data  from  the  literature.  A  new  method  is  proposed  to 
quantify  the  contribution  of  the  reversible  WGSR  to  CO  production 
by  comparing  the  fluxes  of  CO  at  the  cathode-electrolyte  interface 
and  at  the  cathode  surface.  The  results  show  that  the  reversible 
WGSR  can  produce  or  consume  CO,  depending  on  the  average  rate 
of  WGSR  in  the  cathode.  This  is  different  from  the  common  under¬ 
standing  on  H20/C02  co-electrolysis.  The  electrochemical  model 
will  be  extended  to  2D/3D  model  in  a  subsequent  study. 
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2.  Model  development 

2. 1 .  Working  principles 

The  working  mechanisms  of  an  SOEC  for  co-electrolysis  of  H20 
and  C02  are  shown  in  Fig.  1.  In  the  SOEC,  the  gas  mixture  of  H20, 
C02,  H2,  and  CO  flows  in  the  cathode  channel  while  air  flows  in  the 
anode  channel.  In  the  porous  cathode,  both  H20  and  C02  molecules 
diffuse  through  the  porous  electrode  to  the  triple-phase-boundary 
(TPB)  at  the  cathode-electrolyte  interface,  where  they  are  reduced 
to  H2  and  CO  via  reactions  (1)  and  (2),  respectively. 


H20  +  2e-  - 

>  H2  +  02- 

(1) 

C02  +  2e-^ 

CO  +  O2- 

(2) 

The  oxygen  ions  (O2-)  transport  through  the  dense  electrolyte 
to  TPB  at  the  electrolyte-anode  interface,  where  they  lose  elec¬ 
trons  to  form  oxygen  molecules  (Eq.  (3)).  The  produced  oxygen 
molecules  subsequently  diffuse  through  the  porous  anode  to  the 


anode  surface  and  get  collected. 

202-  -*■  02  +  4e_  (3) 

The  overall  reactions  for  H20  electrolysis  and  C02  electrolysis 
can  be  written  as: 

H2O^H2  +  0.5O2  (4) 

C02  ->  CO  +  0.502  (5) 

In  addition  to  the  above  mentioned  electrochemical  reactions, 
reversible  WGSR  also  occur  in  the  cathode  (Eq.  (6)). 

H20  +  CO  ±5  H2  +  C02  (6) 


In  operation,  the  required  potential  (V)  applied  to  SOEC  can  be 
expressed  as: 


V  —  E  +  Pact, a  +  Pact,c  +  P  ohmic  (7) 

where  E  is  the  equilibrium  potential  (Nernst  potential);  p0hmic  is  the 
ohmic  overpotential;  pact,a  and  pact,c  are  the  activation  overpoten¬ 
tials  at  the  anode  and  cathode,  respectively. 


2.2.  Equilibrium  potentials  including  concentration 
overpotentials 


The  concentration  overpotentials  are  not  explicitly  expressed  in 
Eq.  (7)  as  they  are  implicitly  included  in  the  Nernst  potentials  (Eqs. 
(8)  and  (9))  for  reactions  (4)  and  (5),  respectively  [29,46]: 


£H2=E°2  +  gln 


V(po2. 

ph2o 


A/2 


(8) 


£co=E£0  +  §ln 


pl  .(pi  'j1/2 
1  CO  V  02  > 

pl 

C02 


(9) 


where  E°  is  the  voltage  under  standard  conditions;  P^yP^Q,  P!c q2> 
P£0  and  P^2  are  the  partial  pressures  of  H2,  H20,  C02,  CO  and  02 
at  the  electrolyte-electrode  interface,  respectively.  T  is  tempera¬ 
ture  (I<).  R  is  the  universal  gas  constant  (8.3145  J  mol-1 I<-1 );  and  F 
is  the  Faraday  constant  (96,485  C  mol-1 ).  The  E°  can  be  calculated 
from  thermodynamics  ( AG/(2F)).  At  600  K  and  1 200  K,  the  values  of 
E^  are  1.109017  V and  0.940172  V,  respectively  [47].  Similarly,  the 
values  of  £°0  are  1.195502  V  and  0.923869  V,  respectively.  Assum¬ 
ing  linear  variation  of  E°  between  600  K  and  1200K,  the  Nernst 


potentials  (including  concentration  overpotentials)  for  Eqs.  (4)  and 
(5)  can  be  written  as: 


RT 

EEo  =  1.253-0.00024516  T+  —  In 
z  2  F 


RT 

Eco  =  1.46713  -  0.0004527T+  —  In 

2  F 


PH2(V 


Pl 


h2o 


pco(po2) 

pco2 


0.51 


(10) 


(11) 


2.3.  Multi-component  mass  transfer  in  porous  electrodes 


In  order  to  determine  the  gas  partial  pressures  at  the 
electrolyte-electrode  interface,  the  multi-component  reactive- 
transport  processes  in  the  porous  electrodes  must  be  solved 
(Eq.  (12)).  For  modeling  of  multi-component  reactive-transport  in 
porous  media,  Fields  Model  (FM),  Stefan-Maxwell  Model  (SMM), 
and  the  Dusty-Gas  Model  (DGM)  have  been  widely  used  in  the  lit¬ 
erature  [48].  The  DGM  is  used  in  the  present  study  due  to  its  better 
accuracy  for  multi-component  gas  transport  [49].  As  the  reversible 
WGSR  does  not  change  the  total  molar  number  of  gas  species,  the 
pressure  effect  on  mass  transfer  can  be  safely  neglected  [48,50]. 
According  to  DGM,  the  reactive-transport  of  gas  species  i  can  be 
expressed  as  [49]: 


£  a(yfP) 
RT  dt 


=  —VNj  +  Rj 


(12) 


Ni  yjNi  -yjNj  P  dyj 

neff  +  2^  neff  RT  dx 

Uk  j=\,j±i  Uij 


where  yz  is  the  molar  fraction  of  species  i;  s  is  electrode  porosity, 
Rj  is  the  reaction  rate  (molm-3  s-1);  is  the  effective  Knudsen 

diffusion  coefficient  (m2  s-1 )  of  species  i;  is  the  effective  binary 

diffusion  coefficient  (m2  s-1 )  of  species  i  and  j;  Pis  the  pressure  (Pa). 
Nf  is  the  flux  of  species  i  (molm-2  s-1).  dc  and  da  are  thicknesses 
of  cathode  and  anode,  respectively,  x  is  the  depth  inside  the  elec¬ 
trode,  measured  from  the  electrode  surface,  as  can  be  seen  from 
the  computational  domain  shown  in  Fig.  1. 

The  effective  binary  diffusion  coefficient  (D/f )  can  be  evaluated 
as: 


Deff  = 


£  0.0026 T1-5 
!  p^/M^crfj^D 


(14) 


Mij  = 


2 

(1/M;)  +  (1/Mj) 


(15) 


where  s/§  is  the  ratio  of  porosity  to  tortuosity  of  porous  electrodes; 
cfij  is  the  mean  characteristic  length  of  species  i  and  j;  £2D  is  a 
dimensionless  diffusion  collision  integral,  which  can  be  calculated 
as: 


(16) 


£2d 


1.06036  0.193  1.03587  1.76474 

ro.i56i  +  exp(0.47635r)  +  exp(1.52996r)  +  3.89411r 


(17) 


(18) 


Here  kb  =  1.38066  x  10-23(J  K-1 )  is  Boltzmann’s  constant.  The  val¬ 
ues  of  <7,  and  £u  can  be  used  in  the  present  study  are  summarized 
in  Table  1  [51]. 
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h2o,  CO  2  H2,  CO 


X=0 


Cathode 


x  =dt 

Electrolyte 


x=d, 


Anode 


x=0 


•  Ni 
Q  YSZ 
Q  LSM 


co2+h2^±co+h2o 


H20+2e  -±H2+(?- 
C02  +2e~  —>  CO +02~ 


■201'  -+Q  +4e~ 


Fig.  1.  Working  mechanisms  of  SOEC  for  co-electrolysis  of  H20  and  C02. 


The  Knudsen  diffusion  coefficient  can  be  calculated  as: 


Deff  = 


s2rp  [s RT 


(19) 


where  rp  is  the  radius  of  pores. 

In  the  porous  cathode,  the  transport  of  gas  species  is  related 
to  the  rate  of  reversible  WGSR  and  the  rate  of  electrolysis  reaction 
(current  density).  The  rate  ( Rwgsr ,  mol  nrr3  s-1 )  of  reversible  WGSR 
can  be  determined  by  a  widely  used  formula  as  [52-59]: 

Rwgsr  =  Kf  ^Ph2oPco  -  ^H^c°2  ^  (20) 

ksf  =  0.0171  exp  (~1^191)  (mol  m“3  Pa“2  s"1)  (21) 


Kps  =  exp(-0.2935Z3  +  0.6351Z2  +4.1788Z  +  0.3169) 


Z  = 


1000 

W) 


(22) 

(23) 


The  local  mass  conservation  can  be  applied  as: 


dNH2 

dx 


=  Rwgsr I 


dNH2o 

dx 


=  -Rwgsr 


(24) 


The  above  governing  equations  for  describing  the  multi- 
component  mass  transfer  are  inter-related  differential  equations 
and  can  be  solved  by  numerical  method.  The  fourth  order 
Runge-Kutta  method  is  adopted  for  solving  the  DGM  equations. 
After  obtaining  the  partial  pressure  of  gas  species  at  the  TPB 
(electrolyte-electrode  interface),  the  Nernst  potentials  (including 
concentration  overpotentials)  can  be  calculated  using  Eqs.  (10)  and 
(11). 


2.4.  Activation  overpotential 


The  activation  overpotentials  reflect  the  electrochemical  activ¬ 
ity  of  the  electrodes.  In  the  literature,  the  Butler-Volmer  (BV) 
equation  is  the  most  widely  adopted  formula  for  describing  the  acti¬ 
vation  overpotentials  of  SOEC/SOFC.  However,  experimental  works 
suggest  that  the  activation  overpotentials  of  SOEC/SOFC  almost  lin¬ 
early  vary  with  the  current  density  [60].  In  the  present  study,  the 
linear  formula  is  used  [61  ]. 


Pact,  H2,i 


RTJh2 

"h2FJS2,,- 


(29) 


dNco2 

dx 


=  Rwgsr \ 


dNc  o 
dx 


=  -Rwgsr 


(25) 


At  the  cathode-electrolyte  interface,  electrochemical  reactions 
take  place  and  thus  the  flux  of  H2  and  CO  can  be  related  to  the 
current  densities  as: 


Nh2 


I  x=dc 


NCo 


\x=dc 


Jh2  . 

2F  ’ 

Jco . 
2F  ’ 


NH20|x=dc 
NC°2  \X=dc 


Jh2 
2  F 

Jco 
2  F 


(26) 

(27) 


At  the  anode  side,  air  is  usually  used  as  a  sweep  gas.  The  oxygen 
molecules  produced  at  the  TPB  of  anode  transport  to  the  anode 
surface  and  get  collected.  The  flux  of  02  is  related  to  the  current 
density  as: 


Nq2 


\X=da 


Jco  +  Jh2 

4F 


(28) 


Table  1 

Parameters  used  in  calculating  the  effective  diffusion  coefficients  [51  ]. 


CO 

co2 

h2 

o2 

n2 

H20 

CTi 

3.69 

3.941 

2.827 

3.467 

3.798 

2.641 

Silk 

91.7 

195.2 

59.7 

106.7 

71.4 

809.1 

Pact,  CO,  i 


RTJco 


ncoFfc 


CO  ,i 


(30) 


where ] ^  .  and  .  are  the  exchange  current  density  for  H20  elec¬ 
trolysis  and  C02  electrolysis,  respectively.  The  subscript  i  (i  =  a  and 
c)  represents  the  anode  and  cathode,  respectively. 

The  exchange  current  densities  (J^  .  and  f)  represent  the 
readiness  of  the  electrode  to  proceed  with  the  electrochemical  reac¬ 
tion  and  depend  on  the  operating  temperature  as  [62,63]: 

J°2jC  =  k H2,cexp(-^)  ;  Jj  =  kaex P  (~%^)  OD 

Here  /<h2,c  and  ka  are  the  pre-exponential  factors  for  the  cathode 
(for  H20  electrolysis)  and  anode,  respectively.  The  values  of  Eact,c 
and  Eactta  are  1.0  x  105Jmol-1  and  1.2  x  105Jmol-1,  respectively 
[62,63].  Chan  and  Xia  [64]  recommend  the  values  of c  and  at 
1 073  K  to  be  5300  A m-2  and  2000  A m-2,  respectively.  Using  these 
values,  the  pre-exponential  factors  can  be  calculated.  From  exper¬ 
iments,  the  rate  of  electrochemical  oxidation  of  H2  is  found  about 
2.5  times  that  of  CO  [65,66].  Thus,  the  exchange  current  density  of 
the  cathode  for  C02  electrolysis  is  calculated  as  c  =  0.4 /°  c. 
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T  =  1073K 


Current  density,  J  (A.rrf2) 

Fig.  2.  Comparison  of  calculated  data  with  experimental  data  by  Zhan  et  al.  [36]  for 
model  validation. 

2.5.  Ohmic  overpotential 

As  the  interconnector  and  the  electrodes  have  much  higher 
electrical  conductivity  than  the  electrolyte,  only  the  ohmic  over¬ 
potential  of  SOEC  electrolyte  is  considered  in  the  present  study. 
According  to  Ohm’s  law,  the  ohmic  overpotential  of  the  electrolyte 
can  be  calculated  as  [67]: 

Vohmic  =  2.99  X  JL  exp  (32) 

where  L  is  the  thickness  of  the  electrolyte  (m). 

3.  Results  and  discussion 

3.1.  Model  evaluation 

In  this  section,  the  modeling  results  are  compared  with  exper¬ 
imental  data  for  model  validation.  In  the  literature,  Zhan  et  al.’s 
work  [36]  on  syngas  production  by  co-electrolysis  provides  detailed 
experimental  setup  and  operating  conditions,  such  as  the  gas 
composition  at  the  inlet  and  outlet,  the  thickness  of  the  SOEC  com¬ 
ponents,  porosity  of  the  electrodes,  etc.  In  their  study,  the  current 
density-voltage  (J-V)  characteristics  of  C02/H20  co-electrolysis 
by  a  cathode-supported  planar  SOEC  are  measured.  The  thick¬ 
ness  of  Ni-YSZ  (yttria-stabilized  zirconia)  cathode,  LSCF-GDC  (GDC: 
Ceo.9Gdo.1O1. 95 1  LSCF:  Lao.6Sro.4Coo.2Feo.8O3)  anode,  and  YSZ  elec¬ 
trolyte  are  600  p,m,  20-30  [Jim,  and  10[xm,  respectively.  The 
porosities  of  the  cathode  and  anode  are  about  40%  and  30%,  respec¬ 
tively.  At  an  operating  temperature  of  1073  K,  the  molar  fractions 
at  the  inlet  of  the  cathode  are:  50%  H20,  25%  H2,  and  25%  C02.  The 
current  density  and  exhaust  gas  composition  are  measured  at  oper¬ 
ating  potentials  from  0.95  V  to  1.3  V.  These  operating  conditions 
and  structural  parameters  from  the  experiments  [36]  are  used  as 
input  parameters  in  the  theoretical  simulation.  Gas  composition  at 
the  inlet,  exhaust  and  the  average  of  these  two  (should  be  more 
suitable)  are  used  as  gas  composition  at  the  cathode  surface. 

Fig.  2  compares  the  modeling  results  with  the  experimental 
data  from  [36].  The  difference  between  the  modeling  results  and 
experimental  data  is  smaller  when  the  average  gas  composition 
is  adopted.  Although  the  calculated  current  density  is  about  25% 
smaller  than  the  measured  current  density  at  1.0  V,  the  calculated 
data  at  other  potentials  agree  well  with  the  experimental  results. 
In  the  subsequent  parametric  simulation,  typical  data  from  the 


Table  2 

Parameters  used  in  the  simulation. 


Parameter 

Value 

Operating  temperature,  T  (K) 

1073 

Operating  pressure,  P  (atm) 

1.0 

Electrode  porosity,  e 

0.4 

Electrode  tortuosity,  £ 

2.0 

Average  pore  radius,  rp  (pm) 

Cathode-supported  SOEC 

0.5 

Anode  thickness  da  (pm) 

50 

Electrolyte  thickness,  I  (|Jim) 

50 

Cathode  thickness,  dc  (pm) 

500 

Anode  inlet  gas  molar  ratio:  02/N2 

0.21/0.79 

Cathode  inlet  gas  molar  ratio3:  H20/C02/H2/C0 

0.497/0.25/0.25/0.03 

SOEC  operating  potential  (V) 

1.3 

a  Various  gas  compositions  are  studied  and  the  details  can  be  found  from  the  text. 


literature  (i.e.  thickness  of  SOEC  components)  are  used,  as  shown  in 
Table  2.  It  should  be  noted  that  there  are  4  gas  species  in  the  cath¬ 
ode,  although  molar  fractions  of  only  3  gas  species  are  provided. 

3.2.  Effect  of  operating  temperature 

The  effects  of  operating  temperature  on  co-electrolysis  of  C02 
and  H20  in  an  SOEC  are  shown  in  Fig.  3.  It  is  found  that  the  rate 
of  reversible  WGSR  varies  from  negative  values  near  the  cathode 
surface  to  positive  values  near  the  cathode-electrolyte  interface 
(Fig.  3a).  The  results  also  indicate  that  the  WGSR  does  not  reach 
equilibrium.  The  negative  rate  of  WGSR  at  the  cathode  surface  is 
mainly  due  to  relatively  high  molar  fractions  of  C02  and  H2  near 
the  cathode  surface  (Fig.  3b  and  c).  The  positive  rate  of  WGSR  inside 
the  cathode  (particularly  near  the  cathode-electrolyte  interface)  is 
mainly  caused  by  relatively  high  concentration  of  CO  inside  the 
cathode  (Fig.  3c).  It  is  also  found  that  the  rate  of  WGSR  increases 
with  increasing  temperature,  due  to  the  considerable  increase  in 
the  value  of  kSf  from  1.14516  x  10-8  molm-3  Pa-2  s-1  at  873  K  to 
1 .620625  x  1 0-7  mol  nrr3  Pa-2  s-1  at  1 073  K.  The  calculated  rate  of 
reversible  WGSR  is  generally  higher  than  the  data  for  SOFC  with 
internal  reforming,  such  as  [68].  This  is  because  the  rate  of  WGSR 
depends  on  not  only  the  temperature,  but  also  the  products  ofpH2o 
and  pc o  as  well  as  Ph2  and  Pco2>  as  can  be  seen  from  Eq.  (19).  The 
molar  fractions  of  C02  and  CO  in  an  SOFC  fueled  by  pre-reformed 
syngas  are  typically  about  4.4%  and  2.9%,  respectively  [52,68],  much 
lower  than  the  data  used  in  SOEC  for  co-electrolysis.  The  current 
densities  for  H20  electrolysis  and  C02  electrolysis  increase  with 
increasing  temperature,  as  the  activation  overpotentials  and  the 
ohmic  overpotential  are  smaller  at  higher  temperature  (Fig.  3d). 

A  common  understanding  on  co-electrolysis  of  H20/C02  in  an 
SOEC  is  that  the  WGSR  always  contributes  to  CO  production  [38,40]. 
Based  on  existing  experimental  studies,  it  is  difficult  to  confirm 
the  above  expectation  and  difficult  to  quantify  the  contribution  of 
WGSR  to  CO  production  [45].  It  is  understood  that  mass  transfer 
becomes  limiting  only  at  very  high  current  density  and  thus  the 
CO  production  is  governed  by  electrochemical  reaction  (electrol¬ 
ysis)  and  chemical  reaction  (WGSR).  In  this  study,  a  new  method 
is  proposed  to  examine  the  contributions  of  reversible  WGSR  and 
electrolysis  to  CO  production,  by  comparing  the  fluxes  of  CO  at  the 
cathode-electrolyte  interface  and  at  the  cathode  surface  respec¬ 
tively.  The  flux  of  CO  caused  by  electrochemical  reaction  can  be 
easily  calculated  by  Eq.  (27).  The  difference  in  CO  fluxes  at  the  cath¬ 
ode  surface  and  at  the  cathode-electrolyte  interface  is  caused  by  the 
reversible  WGSR.  At  a  temperature  of  873  K,  the  fluxes  of  CO  and  H2 
are  negative  at  the  cathode  surface.  This  means  these  species  diffuse 
from  the  inside  of  the  cathode  to  the  cathode  surface  (Fig.  3e).  Larger 
CO  flux  at  the  cathode  surface  than  at  the  cathode-electrolyte  inter¬ 
face  (Eq.  (27))  means  that  the  reversible  WGSR  contributes  to  CO 
production.  The  reversible  WGSR  and  C02  electrolysis  contribute  to 
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S:  Flux  at  the  cathode  surface 


Fig.  3.  Effect  of  temperature  on  co-electrolysis  of  H20  and  C02  in  an  SOEC  at  an  operating  potential  of  1.3  V  and  inlet  gas  composition  (molar  fraction)  of  H20:  49.7%;  H2: 
25%;  C02 :  (25%);  and  CO  (0.3%)  -  (a)  rate  of  reversible  WGSR;  (b)  molar  fraction  of  H20  and  H2 ;  (c)  molar  fraction  of  C02  and  CO;  (d)  current  density;  (e)  flux  (mol  rrr2  s_1 ) 
of  CO  and  H2  at  873  K;  (f)  flux  (mol  irr2  s-^ofCO  and  H2  at  1073  K. 


CO  production  by  25%  and  75%,  respectively  (Fig.  3e).  Accordingly, 
the  reversed  WGSR  consumes  H2,  leading  to  lower  H2  flux  at  the 
cathode  surface  than  at  the  cathode-electrolyte  interface  (Fig.  3e). 
At  1073  K,  the  CO  flux  at  the  cathode  surface  is  about  19%  lower 
than  at  the  cathode-electrolyte  interface  (Fig.  3f),  meaning  that  CO 
is  consumed  by  WGSR.  The  above  analyses  show  that  the  reversible 
WGSR  can  produce  or  consume  CO,  depending  on  the  average  rate 
of  WGSR.  It  is  also  different  from  the  common  understanding  on 
co-electrolysis  that  WGSR  always  contribute  to  CO  production. 

3.3.  Effect  of  gas  composition  at  the  cathode  surface 

The  inlet  gas  composition  is  varied  to  examine  its  effect  on 
co-electrolysis  behavior.  Three  cases  are  examined  with  different 
molar  fractions  at  the  cathode  surface  -  case  1:  H20  (49.7%),  H2 
(25%),  C02  (25%),  CO  (0.3%);  case  2:  H20  (49.7%),  H2  (0.3%),  C02 


(25%),  CO  (25%);  and  case  3:  H20  (25%),  H2  (25%),  C02  (49.7%),  CO 
(0.3%).  Fig.  4a  shows  the  comparison  of  the  calculated  rates  of  WGSR 
inside  the  porous  cathode  for  the  three  cases,  at  operating  tem¬ 
perature  of  1073  K,  potential  of  1.3  V.  It  is  found  that  the  rate  of 
WGSR  is  positive  in  most  of  the  cathode  layer  for  case  2  -  decreas¬ 
ing  from  about  200  mol  m-3  s-1  at  the  cathode  surface  to  about  0  at 
the  cathode-electrolyte  interface  (Fig.  4a).  For  the  other  two  cases 
(1  and  3)  with  very  low  CO  molar  fraction  at  the  cathode  surface, 
the  rate  of  WGSR  increases  from  negative  values  at  the  cathode 
surface  to  slightly  positive  values  at  the  cathode-electrolyte  inter¬ 
face.  The  molar  fractions  of  gas  species  in  the  cathode  for  cases 
2  and  3  are  shown  in  Fig.  4b  and  c.  High  molar  fraction  of  CO  is 
observed  inside  the  cathode,  indicating  the  slow  diffusion  of  CO 
from  the  inside  of  the  cathode  to  the  surface  of  the  cathode.  For 
comparison,  the  light  H2  molecules  can  diffuse  more  easily  thus 
the  molar  fraction  variation  is  smaller.  Despite  the  different  rates 
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Fig.  4.  Effect  of  gas  composition  at  the  cathode  surface,  at  an  operating  temperature  of  1073 1<  and  an  operating  potential  of  1.3  V  -  (a)  rate  of  reversible  WGSR;  (b) 
distributions  of  gas  molar  fraction  at  an  inlet  gas  composition  (molar  fraction)  of  H20:  49.7%;  C02:  25%;  CO:  25%;  and  H2:  0.3%;  (c)  distributions  of  gas  molar  fraction  at 
an  inlet  gas  composition  (molar  fraction)  of  H20:  25%;  C02:  49.7%;  H2:  25%;  and  CO:  0.3%;  (d)  current  density;  (e)  flux  (molrrr2  s_1)  of  CO  and  H2  at  cathode  surface  gas 
composition  of  H20:  49.7%,  C02 :  25%,  H2 :  25%,  CO:  0.3%;  (f)  Flux  (mol  rrr2  s_1 )  of  CO  and  H2  at  cathode  surface  gas  composition  of  H20:  25%,  C02 :  49.7%,  H2 :  25%,  CO:  0.3%. 


of  WGSR,  there  is  no  big  difference  in  the  calculated  current  den¬ 
sity  for  different  gas  compositions  at  the  cathode  surface  (Fig.  4d). 
The  fluxes  of  CO  and  H2  are  shown  in  Fig.  4e  and  f  for  cases  2 
and  3,  respectively.  For  case  2,  the  CO  flux  increases  from  about 
-0.02  mol  m-2  s-1  at  the  cathode-electrolyte  interface  to  about 
+0.02  mol  m-2  s-1  at  the  cathode  surface  (Fig.  4e).  The  negative  flux 
at  the  electrode-electrolyte  interface  (E  in  Fig.  4e)  represents  CO 
production  by  electrochemical  reaction  (electrolysis).  The  positive 
flux  at  the  cathode  surface  means  that  CO  flows  into  the  cath¬ 
ode,  as  WGSR  consumes  CO.  Since  H2  is  produced  from  reversible 
WGSR,  the  flux  of  H2  at  the  cathode  surface  is  higher  than  at  the 
cathode-electrolyte  interface.  For  case  3,  the  CO  flux  is  consider¬ 
ably  higher  at  the  cathode  surface  than  at  the  cathode-electrolyte 
interface  (Fig.  4f),  due  to  the  high  rate  of  the  reversed  WGSR  inside 
the  cathode. 


3.4.  Effect  of  applied  potential 

The  effect  of  applied  potential  on  co-electrolysis  behavior  of 
the  SOEC  is  shown  in  Fig.  5.  It  is  found  that  the  rates  of  WGSR 
vary  from  about  -200 molm-3  s-1  at  the  cathode  surface  to  pos¬ 
itive  values  at  the  cathode-electrolyte  interface  as  the  potential 
is  increased  from  1.0  V  to  1.4  V  (Fig.  5a).  As  expected,  the  current 
densities  increase  with  increasing  operating  potential  (Fig.  5d), 
which  in  turn  increases  the  molar  fraction  of  CO  and  FI2  inside 
the  cathode  (Fig.  5b  and  c).  At  a  relatively  low  operating  poten¬ 
tial  (1.0  V),  the  CO  flux  at  the  cathode  surface  is  about  3  times 
of  the  CO  flux  at  the  cathode-electrolyte  interface  (Fig.  5e),  indi¬ 
cating  the  significant  contribution  of  the  reversed  WGSR  to  CO 
production.  In  addition,  the  positive  H2  flux  at  the  cathode  sur¬ 
face  means  the  FI2  molecules  diffuse  into  the  cathode  and  are 
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e  S:  Flux  at  the  cathode  surface 


Fig.  5.  Effect  of  operating  potential  at  an  operating  temperature  of  1073  K  and  inlet  gas  composition  (molar  fraction)  of  H20:  49.7%,  C02:  25%,  H2:  25%,  CO:  0.3%  -  (a)  rate 
of  reversible  WGSR;  (b)  molar  fraction  of  H20  and  H2;  (c)  molar  fraction  of  C02  and  CO;  (d)  current  density:  (e)  flux  (molirr2  s-1)  of  CO  and  H2  at  an  operating  potential  of 
1.0  V;  (f)  flux  (mol  rrr2  s-1)  of  CO  and  H2  at  an  operating  potential  of  1.4  V. 


consumed  by  the  reversed  WGSR  (Fig.  5e).  At  an  operat¬ 
ing  potential  of  1.4  V,  the  flux  of  CO  at  the  cathode  surface 
is  smaller  than  at  the  cathode-electrolyte  interface  (Fig.  5f), 
indicating  the  average  rate  of  WGSR  is  positive.  From  the 
above  analysis,  it  can  be  seen  that  the  reversible  WGSR  can 
contribute  significantly  to  CO  production  at  a  low  operating 
potential  but  it  can  consume  CO  at  a  high  operating  poten¬ 
tial. 

4.  Conclusions 

An  electrochemical  model  is  developed  to  characterize  the  per¬ 
formance  of  an  SOEC  used  for  H20/C02  co-electrolysis.  The  model 
is  validated  by  comparing  the  simulation  results  with  experimental 
data  from  the  literature. 


It  is  found  that  the  rate  of  WGSR  in  the  porous  cathode  can 
be  positive  or  negative,  depending  on  the  temperature  and  local 
gas  composition.  A  new  method  is  proposed  to  quantify  the  con¬ 
tribution  of  reversible  WGSR  to  CO  production,  by  comparing  the 
CO  fluxes  at  the  cathode-electrolyte  interface  and  at  the  cathode 
surface.  It  is  also  found  that  reversible  WGSR  can  produce  CO  or 
consume  CO,  depending  on  the  average  rate  of  WGSR  in  the  porous 
cathode.  In  the  simulated  case  at  an  operating  temperature  of  873, 
the  reversible  WGSR  contributes  to  about  25%  of  CO  production, 
while  at  1073K,  the  reversible  WGSR  consumes  about  19%  of  CO 
produced  from  C02  electrolysis.  In  addition,  the  gas  composition 
at  the  cathode  surface  greatly  influences  the  distribution  of  WGSR 
rate  inside  the  cathode.  Large  negative  rate  of  WGSR  is  observed 
near  the  cathode  surface  with  very  low  molar  fraction  of  CO.  Gen¬ 
erally,  the  WGSR  plays  an  important  role  in  CO  production  at  a  low 
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operating  potential  while  the  WGSR  can  consume  CO  at  high  oper¬ 
ating  potential  at  an  operating  temperature  of  1073  K  and  inlet  gas 
composition  (molar  fraction)  of  H20:  49.7%,  C02:  25%,  H2:  25%,  CO: 
0.3%. 

The  present  study  provides  detailed  information  for  better 
understanding  the  working  mechanism  of  SOEC  used  for  H20/C02 
co-electrolysis.  The  new  method  proposed  is  useful  to  quantify 
the  contributions  of  WGSR  and  C02  electrolysis  to  CO  production. 
The  electrochemical  model  can  be  integrated  into  2D/3D  model  for 
more  detailed  analysis  and  simulation. 
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